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Abstract. A simplified frequency integrated radiative 
transfer equation is solved to study Compton scatterings 
in the corona of the disk by using numerical iterating 
method. We find that the vertical thickness of the corona 
cannot be used as the typical length to measure the op- 
tical depth of the corona. A semi-analytical approach is 
proposed to calculate the energy dissipations in the corona 
of the disk. We demonstrate that our approach can repro- 
duce the numerical solutions to an accuracy of < 2%. 
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1. Introduction 

' The so-called UV bump and the X-ray power law spectra 
are the most important features in the spectra of Active 
Galactic Nuclei (AGN). ft is well believed that the lumi- 

[ nosities of the AGN are produced by accretion of mat- 
ter through disks onto massive black holes. Shields (1978) 

^ suggested that the UV bump observed in AGN could be 
attributed to the thermal radiation from the disk. Malkan 
(1983) showed that the UV bump in several quasars could 
be fitted with predicted spectra of optically thick accretion 
disks. The observed spectrum in the medium X-ray range 
is close to a power law, with a small dispersion in the val- 
ues of the spectral index, whose average for Seyfert galax- 
ies is 0.7 (Mushotzky 1984; Turner & Pounds 1989). In 
fact, the average X-ray spectral index becomes ^ 0.9 af- 
ter subtracting the component of X-rays reflected by cold 
matter (Pounds et al. 1990). 

The standard geometrically thin, optically thick ac- 
cretion disk model (Shakura & Sunyaev 1973) can be 
employed to explain the UV bump (Malkan 1983), but 
the standard disk is too cold to produce the X-ray radi- 
ation. Another solution for the inner parts of accretion 
disks around black holes is the optically thin, hot disk 
model, which was first proposed by Shapiro, Lightman & 
Eardley (1976). Electron temperature of the disk in this 



model becomes ^ 10^ K and it can explain observed X- 
ray emissions. However, it fails to reproduce the observed 
UV bump. Haardt & Maraschi (1991) proposed a two- 
phase thermal disk-corona model. In this model, most of 
the thermal soft photons emitted from the optically thick, 
cold disk pass the optically thin corona without being scat- 
tered and are observed as UV bump. Only a small fraction 
of them is Compton upscattered to the X-ray range by the 
high temperature thermal electrons in the corona. About 
half of the Compton scattered photons are directed to- 
ward the disk and are reprocessed to emerge as blackbody 
radiation, while the remaining half are directed upwards 
and are observed as power law X-ray spectrum. In this 
model, the gravitational energy of the accreting matter is 
mainly released in the corona of the disk. Along the line 
of this model, some workers investigated the structure of 
the disk-corona system with different assumptions (Naku- 
mura & Osaki 1993; Kusunose & Mineshige 1994; Svens- 
son & Zdziarski 1994). The X-ray spectra emitted by the 
corona are studied in detail by using different approaches, 
especially the Monte Carlo methods (Haardt 1993; Haardt 
& Maraschi 1993; Titarchuk 1994; Poutanen & Svensson 
1996). Nakamura & Osaki (1993) found in their numeri- 
cal simulations that the effective optical depth of photons 
traveling through the corona may be larger than the op- 
tical depth in the normal direction of the disk. 

In this work, we solve the integro-differential radiative 
transfer equation by using numerical iterating method. 
The basic equations describing the problem are listed in 
Sect. 2. Section 3 contains the results and a semi-analytical 
approach proposed to reproduce the numerical results. 
The last section is the discussion of the results. 

2. Basic equations 

The frequency-integrated radiative transfer equation is 

a/(T,/i) 



dr 



/(r,/i)-5(T), 



(1) 
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where r is the optical depth measured in the normal di- 
rection of the disk, /(t, /i) and S{t) are the frequency- 
integrated specific intensity and source function, fj, = 



2 
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cos 9, describes the direction of 7(r, /u) with respect to the 
normal direction. 

Here we only consider pure Compton scatterings in 

the corona and assume the Compton scatterings to be 
isotropic. Thus, the source function S{t) could be written 
as 



s{t) = ^^ I i{T,ij)dQ. = J i{T,fx)d,,, 



(2) 



where A = 40 -|- 160^, is the mean fraction of the energy 
change of the photon through single scattering. The di- 
mensionless electron temperature 6 = kTe/meC^, where 
k, me and c are the Boltzman constant, rest mass of the 
electron and the light speed. Strictly speaking, the Comp- 
ton scattering in the corona is anisotropic (Ghisellini et 
al. 1991; Haardt 1993). Here we only intend to investigate 
the total energy transferred from the mild relativistic elec- 
trons in the corona to the scattered photons. In the first 
attempt, we employ Eq. (2) to describe the source function 
in the radiative transfer equation for the sake of simplic- 
ity. We believe this will not affect the main results of our 
present investigation. 

We consider the homogeneous isothermal layer in a 
plane-parallel geometry, which can model the hot corona 
over an accretion disk. Compton parameter y is defined 
by 



y = To(4e + 16e2)max(l,ro), 



(3) 



where tq is the optical depth of the corona by electron 
scattering measured in the vertical direction. Now, we can 
solve Eq. (1) combining Eq. (2) with appropriate bound- 
ary conditions, if the electron temperature Q is specified. 

In this work, we assume the input soft photons from 
the cold disk are isotropic, then the boundary conditions 
for the problem are as follows: 



J(r, /x) = Jo, M > 0, at t = tq 
and 

/(r, = 0, At < 0, at r = 0, 



(4) 



(5) 



where t is the electron scattering optical depth measured 
in the normal direction of the disk. The energy flux of the 
soft photons from the cold disk is 



Fs = 2tt j I{To,ii)ndiJb = tt/q. 



(6) 



Most of the soft photons emitted by the cold disk pass the 
optically thin corona without being scattered. The flux of 
the soft photons escaped from the system is given by 



1 

Fesc = 2Ty j I{To,fJ,)e-^°/^f^dfi. 



(7) 



A small fraction of soft photons from the cold disk is 
Compton upscattered in the corona. About half of the 
scattered photons are directed upwards. The remaining 
half are directed downwards, part of them are absorbed 
by the cold disk and part are reflected. 

The fluxes of the upward scattered photons and down- 
ward scattered photons are 



Fuc = 27r 



and 



Fdc = 27r 



± 

j I{0,iJ,)iJ.diJ.- 



V 

j I{To,n)iJ,diJ,, 



(8) 



(9) 



respectively. 

The fraction D of the energy dissipated in the hot 
corona to the input soft photon energy flux from the cold 
disk is then given by 



D 



FuC + Fac + F, 



esc ^ 8 



(10) 



Finally, the fraction t] of Compton upscattered photon en- 
ergy directed downward to the cold disk is 



FdC 



Fdc + Fuc 



(11) 



3. Results 



Equation (1) is an integro-differential equation and cannot 
be solved by simply iterating, since the source function is 
not known a priori and depends on the solution 7(t, jj) at 
all directions through a given point. In our mimcrical cal- 
culations, we first set the source function S{t) = and the 
test solution /i(r, /x) is available. Let /(r, /x) = /i(t, /x), 
the source function 5*1 is obtained from Eq. (2). Then 
the more accurate solution I2 (r, fi) could be obtained by 
iterating Eq. (1), if we let S{t) = 5i(r) in Eq. (1). Sim- 
ilar to the first step we have taken, the solution /„(t, /i) 
and the source function S'„(r) will be known by repeat- 
edly iteration. The iteration is terminated until it reaches 
{Sn - Sn-i)/Sn < 10~^, then the solution to Eq. (1) is 
available. 

The numerical results on the fraction D of the energy 

dissipated in the hot corona to the input soft photon en- 
ergy fiux from the cold disk are plotted in Fig. 1. It is found 
that the behaviour of D varies with the different electron 
temperatures Q. For the higher electron temperature, the 
value D becomes larger even at the same Compton pa- 
rameter y. 

Considering that the soft photons only suffer single 
scattering in the corona, we can obtain the mean energy 
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Fig. 1. The numerical results of the fraction D of the en- 
ergy dissipated in the corona to the input soft photon 
energy vs. Compton parameter y, with respect to the elec- 
tron temperatures = 0.1, 0.2, 0.3, 0.4. 

gained by the photons from the hot electrons in the opti- 
cally thin corona is 

1 

AE = A [ /(to,m)(1 - e--"/^')fidi^ = loy^, (12) 

J To 


where 

1 

a(ro) = /(I - e-^°^ni^dfi, (13) 



is the fraction of the photons which suffer single scattering 
in the corona to the injecting isotropic soft photons. It can 
be seen in Fig. 2 that the value a(To)/To decreases with the 
tq. Equation (12) indicates that the mean energy gained 
by soft photons through first scattering depends not only 
on the Compton parameter y, but also on the optical depth 
To of the corona. Only when y <^ 1, the mean energy 
gained per scattering could be written as oc y. It is shown 
from Eq. (3) that the optical depth becomes smaller for 
the higher electron temperature of the corona when the 
value of y is fixed. Prom Eqs. (12), (13) and Fig. 2, we 
know that, for the same Compton parameter y, there is 
more energy gained by the soft photons through single 
scattering for the corona with a smaller optical depth tq, 
i.e., a higher electron temperature Q. This is the reason 
why the D — y behaviour depicted in Fig. 1 varies with 
the electron temperature of the corona. 

The relation _D '--^ 2 is required in the disk-corona sys- 
tem (Haardt & Maraschi 1991; Nakamura & Osaki 1993). 



Fig. 2. The function a(ro)/To vs. vertical optical depth tq 
of the corona. 

The Compton parameter y is about 0.6 (Haardt & Matt 

1993). Our numerical results show that y 0.35 — 0.4 is 
more desirable to yield D ^ 2 corresponding to different 
electron temperatures O. The difference may be attributed 
to an ambiguity in measuring the optical depth of the sys- 
tem with slab geometrical configuration. For the photon 
travels just in the normal direction of the slab and the 
optical depth tq <C 1, the Compton parameter y well de- 
scribes the situation. Nonetheless, most photons pass the 
slab in the directions other than the vertical direction, 
these photons will travel a longer distance than the ver- 
tical length before they leave the corona since the disk is 
in slab geometry. The effective optical depth of these pho- 
tons arc larger than tq. Hence the Compton parameter 
y defined by Eq. (3) cannot describe the corona in slab 
configuration correctly. 

It should be indicated that the present calculations 
are not effective in the case that saturation of Compton 
scatterings is important. The typical blackbody radiation 
temperature of the cold disk in Active Galactic Nuclei is 
~ 5 — 50 eV (Haardt & Maraschi 1993). Assuming the 
temperature of electron in the corona 6 ~ 0.5, the mean 
scattering number of the photon is 5 before the Comp- 
ton process saturates. For a lower electron temperature O, 
the photon suffers more scatterings before it reaches sat- 
uration, for example, the mean scattering number is > 20 
for 6 ~ 0.1. Hence, the saturation of the Compton pro- 
cess in the disk should be taken into account at least in 
some cases with relatively higher electron temperatures. 
In principle, the recoil of the electron in the scattering 
and Klein-Nishina electron scattering cross-section should 
be taken into account in the study of the saturation of 
the Compton process for the mild relativistic electrons 
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treated here. For simplicity, we do not apply the com- 
plicated formalisms including Klein-Nishina cross-section, 
instead, we just repeat the calculations in the same way 
described at the beginning of this section till the mean 
energy of the scattered photons exceeds that of the elec- 
trons in the corona, and then let it simply be the mean 
energy of the electrons. Yet, this numerical approach is 
still too complicated to be used in constructing the disk- 
corona model. So, we propose a semi-analytical approach 
to approximate the numerical results. 

Suppose the number of photons which suffer at least 
fc-fold scatterings is Nk, the mean energy gained by the 
photons in the fc-th scattering is Nk{l + A)''^^Aei, where 
ej is the mean energy of the input soft photons from the 
cold disk, A = 4G -|- 160^, is the amplified factor. Thus, 
the fraction D of the energy dissipated in the hot corona 
to the input soft photon energy flux from the cold disk is 



(14) 



fe=i 



where m is the scattering number of the photon before 
the Compton process saturates, A^o is the number of in- 
put soft photons from the cold disk. We know that the 
effective optical depth will be larger than the vertical op- 
tical depth tq. The fraction of the soft photons suffering 
the first scattering is a(To) given by Eq. (13). The first 
scattered photon number Ni = Noa{To). We further as- 
sume Nk+i = Nkb{To) for k > 1, where the coefficient 
6(to) is assumed to be a constant for the multiple scat- 
terings except the first one for the photons. Then, 6(to) 
is only a function of the optical depth tq. Therefore, Eq. 
(14) could be rewritten as 



D = 



a{To)A[l-b{To)"'{l + Ay 

l-b{To){l+A) 



(15) 



We take Eq. (15) as the 'seed' analytic form to fit the 
numerical solutions and find that 



D = 



a{To)A[l~b{Toni + Ar] 
1-6(to)(1+A) 



(16) 



can well approximate the numerical results to an accuracy 
of < 2% in the range of Z? = — 5 for the optically thin, 
hot corona. The two coefficients a(To) and 6(ro) in the 
formula are given by Eq. (13) and 



6(ro) = (0.2302ro°'^'^^^ + 0.5102ro) 



a(To) 
To 



(17) 



respectively. The mean eff'ective scattering number x for 
the photon before it reaches saturation is 



log(i 



0.511 X 10*^6 
kJlb 



0.: 



(18) 



where 0.8 in the right side of the equation is simply in- 
duced to fit the numerical results better, kTi,), (in the 



Fig. 3. Comparisons between the results of the semi- 
analytical approximations and the numerical results 
(small circles) for the blackbody radiation temperatures 
of the soft photons kT^i, = 0, 5 and 50 eV, with respect to 
the electron temperature 9 = 0.15. 



Fig. 4. Same as Fig. 3, but the electron temperature G 
0.45 instead. 



unit of eV) is the temperature of the blackbody radiation 
from the cold disk. The comparisons between the semi- 
analytical approximations and the numerical results of the 
problem are plotted in Figs. 3 and 4. 
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4. Discussion 

Wc investigate the energy dissipations in the hot corona 
of the disk by iterating the radiative transfer equation. A 
somewhat simpUfied method is employed to study the sat- 
urations of Compton processes. In Fig. 3 both the numeri- 
cal results and semi-analytical approximations are plotted. 
We find that the value of D decreases with the energy of 
the soft photons kT^i, for the same electron temperature 
Q. The reason is that in the limit of y ^ 1, all photons are 
in equilibrium with the electrons, D ^ 0.511 x lO^S/fcT^b. 
Therefore, a larger fcTbt, for the input soft photons im- 
plies a smaller D for the fixed electron temperature Q. 
We note that the behaviours oi D — y do not change sig- 
nificantly with the soft photons fcTf,;, for relatively lower 
electron temperatures (see Fig. 3 for G = 0.15). In this 
case, the photon suffers more scatterings before it reaches 
satiiration, which indicates only a smaller fraction of in- 
put soft photons is finally in energy equilibrium with the 
electrons in the corona. The D ^ 2 is desired in two-phase 
model, which requires that the Compton parameter y to 
be 0.3 — 0.4 with respect to the different electron temper- 
atures © and the blackbody radiation temperatures kT^b 
of the photons from the cold disk. For relatively higher 
electron temperature, there is a larger amplified factor A, 
and the Compton process will soon saturate within sev- 
eral scatterings except fcTb;, ^ 0. Compared with the lower 
electron temperature case, more soft photons are driven 
to equilibrium with the electrons, and then no longer gain 
energy from hot electrons through further scatterings. We 
suggest that Eq. (16) combining Eqs. (13), (17) and (18) 
could be used as an approximate approach in the study 
of the structure and the spectrum of the disk-corona sys- 
tems. Finally, we point that the pair production and an- 
nihilation in the hot corona should be taken into account 
when the electron temperature in the corona exceeds 10^ 
K. This will be the scope of the further work. 
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